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ABSTRACT 

We discuss decay of unstable particles and pair annihilation of stable heavy parti- 
cles that occur in the cosmic medium, from the point of the fundamental microscopic 
theory. A fully quantum mechanical treatment shows that the effect of thermal en- 
vironment on these processes cannot be described in terms of quantities on the mass 
shell alone, thus requiring an extention of the Boltzmann-like equation. The off-shell 
effect tends to prolong physical processes that take place subsequent to the decay. 
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1. Introduction 

Any important event in cosmology is a series of elementary processes that occur 
in the expanding universe. Some of these events can be understood in terms of 
a rather simple process or their combination. One should however keep in mind 
that the process does not occur in isolation, unlike a well-prepared experiment in 
laboratory. In the past the effect of environment on elementary processes has been 
dealt with in very simplifed manners in practical computation in cosmology, and in 
many cases has been completely ignored. 

We address this problem here and would like to fully incorporate the quantum 
mechanical principles. A subtle quantum mechanical effect might require a change 
of the cherished practice in this field, but we would like to think of the problem, 
going back to the basic point. It turns out that the familiar Boltzmann equation 
must be reexamined due to its neglect of the off-shell effect. PJ 

Although no specific application is discussed in any detail here, it is presumably 
helpful to mention at the outset what I have in mind for applications. They are 
X-boson decay for baryogenesis, and WIMP or LSP annihilation for dark matter. 
In all these cases it is important that decay or annihilation products make up a 
part of the cosmic environment. Thus, when heavy particles are about to decay 
or decouple, inverse processes must be considered simultaneously. Interaction with 
thermal environment becomes the major concern for this problem. 

A usual computational practice of the rate calculation for a particular pro- 
cess goes something like this; one first computes the elementary S-matrix element, 
(f\S\i) , between the initial state % and the final state / on the mass shell, and 
then averages the probability in a thermal bath, ((/IS)?)! 2 . This rate is used in the 
Boltzmann equation in the expanding universe. 

I can think of two basic problems in this approach; first, the states chosen \i) and 
\f) may not be necessarily on the mass shell, as will be explained more fully below, 
and secondly, even the on-shell matrix element 



may contain propagation effects in thermal environment, different from those in 
vacuum. 

I shall first repeat the argument already made on the first point. || For this it is 
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useful to look at the thermal Green function, 

G*(u) = A(u + ie(u))-iJ^, (1) 
with T — 1/(5 the temperature and the spectral function is given by 



(uj 2 -uo 2 k - mi{u)) 2 + ( $m(u) ) : 



= — —2 «m/..^2 i fC^TJf. .U2 > ( 2 ) 



with ujk = Vk 2 + M 2 . Here A is the Green function in vacuum. The proper self- 
energy II(u;) is in general a complicated object, and only in the weak coupling limit 
one may use the Breit-Wigner function for p(u), leading to the familiar form, 

If one further assumes that the relevant a; integral fully includes the peak region 
~ Uk of width T, then the Boltzmann equation follows. 

The off-shell effect appears in two ways; first, in the effective cutoff of the u region, 
oj > uj c (the threshold and taken to be < Uk) for u integral. If the temperature is 
below the threshold, T < u c , the dominant region that contributes, oo < T, misses 
the peak region. Second, if the coupling is not very weak, one must use the correct 
form and not its Breit-Wigner approximation of the the spectral function p(ou). For 
instance, the behavior of QTI(u;) near the threshold u c may become important at low 
temperatures. A systematic treatment of the off-shell effect overcoming these points 
is called for. 

Let us make clear the fundamental difference between the on-shell and the off- 
shell contributions. In elementary quantum mechanics the off-shellness appears in 
virtual states in perturbation formulas. It deals with the short-time behavior and 
thus represents any transition that does not conserve energy. If one considers a decay 
of unstable system of mass M and decay width T, then the off-shell effect is essential 
at times of t ~ 1/M while the on-shell S- matrix element well describes the long-time 
behavior of t > 1/T (the lifetime). As is well known, the on-shell behavior is given 
by the S-matrix, while the off-shell effect must be discussed in terms of the more 
fundamental quantity, the Green function. 
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2. Model of unstable particle decay 

We consider here the case in which one can clearly separate the small system 
in question and the cosmic environment. This means that one should be able to 
ignore mutual interaction among unstable particles and also to ignore the change 
of the environment due to the interaction with the unstable particle. Thus, all the 
nonlinear effect and the back reaction is ignored. 

Within this linear approximation one can think of a harmonic environment of 
arbitrary spectral distribution as well as a bilinear interaction between the small 
system and the environment. We then take the following model Hamiltonian, 

H = H c + H b + H mt 

= w fc c t c+ / duutf(u)b(u) + j duoJa(uj) (b ] {uo)c + c ] b(uj)\ . (4) 

Here c 1 , c are creation and annihilation operators of the unstable particle of energy 
Uk = Vk 2 + M 2 , and tf(u) , b(u) are those of the environment state, usually taken to 
be two-body states of the decay product. The threshold u c = \jk 2 + {mi + m 2 ) 2 and 
mi is the masses of decay product. The instability condition is imposed; uj c < uj k . 

Although for definiteness we consider the unstable particle decay, one may choose 
with some generalization a pair of heavy stable particles for an initial state. 

It should be reminded of that even in this linear model the irreversibility exists; 
it can occur because an unstable particle couples to an infinitely many degenerate 
environment states. A great virtue of this model is that it can be analytically solved 
and in principle one has many explicit results that can test approximation schemes 
one always needs in practice. We shall discuss extention of this model later. 

A way to make integrability explicit is to present the operator solution. One can 
find variables that diagonalize the Hamiltonian; 

I i x r°° \l cr(oj)a(oj') , \ 

B\uj) = b\u) + F(u + i0 + ) - Jc{u) ct + / du V , \' &V) , (5) 

\ v Jlu c U)' — LU — l0 + I 

where 

/•oo 

H= duouoB^(u)B(u) . 

The transformation from c,b(u) to B(u) (and their conjugates) is canonical. Even 
the operator inversion can be written down. The important function F(z) is analytic 
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except along the cut, z > lu c , and explicitly 

F(z) = pr . (6) 

- z + u k - L°° dw 2^ 

We do not explain here the frequency and the wave function renormalization, for 
which we refer to our paper. 

The essential reason why this model is integrable is saturation of unitarity by 
the elastic one. Analogy to the scattering problem is useful and is indeed made very 
evident by considering the energy eigenstate, 

|w> 5 = B t (o;)|0>, 

where |0) is the perturbative (and also the true) vacuum with c|0) = 0. Clearly, 
this state \u)s is a linear superposition of c^O) and tf{uj)\Qi) . When u>k > uj c , the 
original one-particle state c^O) becomes a resonance, or the unstable particle, which 
is \uj)s- The overlap probability 

\(0\c\uj) s \ 2 = H(uj) 

is very fundamental. It is related to the discontinuity along the cut, 

F(u + i0 + ) - F(u - i0 + ) = 2m a(u)\F(u + i0 + )\ 2 = 2m H(cu) . (7) 

In terms of the spectral function 

H(uj) 



(u-u k + m(u)) 2 + (7ia(u)) 2 ' 
where 



IQJ C UJ Z 

From the discontinuity structure it follows that a simple pole exists in the second 
Riemann sheet extended via this discontinuity formula; the original state pole on 
the real axis at u = u>k moves into the second sheet and acquires an imaginary part, 
Ri T = 2mr(E r ) where E r is the real part of the pole (in the weak coupling limit, 
E r pa ujk after renormalization). This pole describes the exponential decay law; the 
non-decay amplitude pa e - iuJ kt-rt/2 _ 

It can also be shown || that this quantity H(uS) is directly related to the Lehmann 
spectral function pl(u) by 

(uj + M) 2 
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The decay law of the unstable particle state |1) = c^O) in vacuum is precisely 
described as follows. One notes that the non-decay amplitude, 

(l\e- tHt \l) = (0\ce- iHt J\0), 

is calculable if one finds the Heisenberg evolution, c(t) = e lHt c e~ lHt . Since the 
operator relation is known, this is given by 

c(t) = - j°° du ^\J)F(uj + iO + )e~ iu}t B(lo) . (9) 
It is useful to rewrite the Heisenberg operator in terms of the original variables; 

c\t)=g(t)J + i I duuJa(uu)h(uj,t)e lu)t b\uu). (10) 
The basic function here is g(t); 

poo 1 r 

g(t)= dcoH(uj)e tuJt = — dzF(z)e lzt , (11) 

Jw c 2m JCo+Ci 

h(to,t)=iF*(cu + iO + )-ik(uj,t), k(u ,t) = i e~ iuJt g(t) . (12) 



The boundary condition for k (cu , t) is k{uj , t) — > as t — > oo, hence 

I r- e i(z-ui)t 
ICo+d 

The non-decay amplitude is given by 



1 r- „i(z-ui)t 

k(u,t) = / dz F(z). (13) 



(l\e- iHt \l) = g*(t) = / dua(u) \F{uj + i0 + )\ 2 e~ iu)t . (14) 

The second form of the u integral for g(t) along Co + C\ makes separation of the 
pole term explicit; the contour Cq encircles the pole (conjugate to the one already 
discussed, due to e tzt factor instead of e~ lzt ) in the second sheet, while C\ consists of 
two lines parallel to the imaginary axis, one in the first and the other in the second 
sheet. It is clear that the C\ integral gives the non-exponential decay law which 
is important both at early and late times of the decay. For instance, at very late 
times the important part of oo integral comes from the threshold region, and if one 
parametrizes this region using the rate T, 

"(^ hn • (15) 



the non-decay amplitude is given by 

~ 2~o ioir 



(l\e- iHt \l) » — i— ^— 77tt~tt e-™'*-*™ /2 , (16) 
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with Q = uj k — uj c . The decay at late times thus follows the power law; 

~ (_L)3(Qt)-2(«+l) . 

3. Decay in cosmic medium 

Amuzingly, this model and the basic technique outlined is taken over to derive 
the decay law in thermal medium. M, || In the thermal medium one is interested in 
time evolution of the occupation number, f(t), to be given by the thermal average 
over a given mixed state, 

f(t) = (c\t)c(t)) i =ix(c\t)c(t)p i ), (17) 

where p, is the density matrix of the mixed state. 
We take for pi the thermal density matrix, 

Pi = e -/»(ff»+*»)/(tr e -P(H>+Hc)^ 

where Hb and H c are the environment and the subsystem parts of the Hamiltonian. 
It is perhaps necessary to explain why we take for the initial mixed state this one 
and not the one like e~^ u using the total Hamiltonian H . The cosmic environment 
we consider is not a stationary state; it adiabatically changes with the cosmic expan- 
sion. With the momentum redshift a time lag occurs during the decay between the 
unstable particle and its lighter decay products. Light decay products have their own 
stronger interaction not written in the Hamiltonian above, and they tend to keep 
in thermal equilibrium with the well-known temperature variation, T(t)a(t) = con- 
stant, with a(t) the cosmic scale factor. On the other hand, the decay and its inverse 
process cannot keep pace with the cosmic expansion. We model this circumstance by 
postulating that the environment by itself is always in thermal equilibrium, while the 
unstable particle starts to lose its thermal correlation with the environment. This is 
the reason we take the initial density matrix uncorrelated between the unstable and 
the decay product particles, although at much higher temperatures they may be in 
thermal contact. Thus, we take for the initial state average 

(ctc) = /(0), (0+6(0;)), = 0, (tf(u)b(u% = 5(u-u')-J—. (18) 
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For the moment we ignore the effect of the cosmic expansion. Time evolution of 
the occupation number f(t) can be derived using the explicit form of c^(t), eq. (|10|) ; 

m = \ 9 (!)?m+r^ ° {w) }^f ■ (i9) 

Its time derivative can be written in the following form, 

| = -r(t)b(t)| 2 /(o) + £rf^M m ltTP > ( 20 ) 

where the time dependent decay rate is defined by 

r «=- 2s !=-! in i 9 ' 2 - < 2i > 

The rate T(t) is time independent during the pole dominant epoch, but at late times 
it decreases as T(t) — > , using the spectral function having the threshold 

behavior fll5|) . An equivalent form for the time evolution equation is obtained by 
eliminating the initial value (/(0)) dependence; 

f% , , 2^(gh*e-^ t )+T(t)\h\ 2 

This form is more convenient to application to cosmology. 

The first form of time evolution equation (|20|) has a clear interpretation if one 
retains only the on-shell contribution; by using the pole dominance for g(t) = 
e iu k t~rt/2 ^ fo e equation is reduced to 



T t = - re '"{f^-^i)- (23 > 

There are oscillatory terms such as e~ r '/ 2 x (cos or sin(co> — u^t) , which averages 
out in the u integral, leading to this result. Thus, if the initial occupation /(0) is 
chosen to be that of the thermal one, then there is no change, while its deviation 
from the equilibrium value is relaxed back by the rate T. 

The second form of evolution equation fl22|) suggests that irrespcetive of the pole 
dominance the occupation number has an asymptote; 

f°° H(uj) f°° <t(oj) 1 

fco = L T^^l ~ Jk (^-^) 2 + (7raM) 2 ^T' (24) 

since both of g(t) and k(t) in h(u> ,t) asymptotically vanish and 

h{uj,t) -> iF*(to + iO + ) . 
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This asymptotic form reveals the essence of our result; only when the full Breit- 
Wigner region for H(u) contributes, the on-shell result gives a good description for 
the asymptotic occupation number. Otherwise, ^ e $X-\ (on-shell result). 

The particle number density is calculated by summing the occupation number 
over independent momenta; 

n(t) = — / dkk 2 f(t,k). (25) 



2tt 2 

Let us work out its asymptotic expression n^, using for f(t,k). When the pole 
term dominates, this gives a well-known result; for the non-relativistic region, 

Uoq «(^)3/2 e -M/T_ 

V 2tt 

On the other hand, at low temperatures the threshold region dominates in the uj 
integral, and one has a power of temperature; 

n ~A(a)^(^T* ^ fa) _ C(a + 4)r(q + 4)r(f + 1) 

Thus, the low temperature behavior is much enhanced compared to that given by 
the on-shell contribution alone. 

The effect of the cosmic expansion is readily incorpolated for the evolution of the 
number density. The easiest way is to write it for the dimensionless yield Y = 
using the dimensionless time variable, r = Ft ; 



Y 

J^3 



7( t , k) = - 2 » _ m Mt , k ) +91 (t, k ))=- • (29) 



The momentum dependent spectral is a(u,k) = ^- a(y to 2 — k 2 ) , where is 
the one in the rest frame. The usual time-temperature relation T oc 1/ should be 
used in this equation. 

The asymptotic solution for this equation is readily found, ignoring the time 
dilatation effect (M/cuk — > 1); 

y _> (^ + 2)C(a + 4)r(q + 3)r(f + 1) J_ T +1 

87r 2 0Fr(f + |) M l M J ' 1 J 
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The important point is that it has a power behavior of the cosmic temperature, as 
is guessed from n M of eq. (|26|) . The power a + 1 is related to the threshold behavior 
of the spectral function a{u). 

I refer to our recent paper || on the detailed time evolution, but roughly Y follows 
initially the thermal behavior, (^pf) 3 ^ 2 e~ M l T , and after a short transient epoch it 
approaches the asymptotic form above. The transition epoch is characterized by two 
temperature scales, T* and T eq given by 



/ V ~ = / 45 m pl F /Ol\ 

M ~ y(a + 4) lnf ' V \ W N M 2 ' 1 ' 



!f»^rj , (32) 

where iV is the number of massless particle species contributing to the energy density. 
The quantity rj is the decay rate relative the Hubble rate calculated at T = M. Usu- 
ally, T* > T eq . The higher temperature T* is given by the time when the exponential 
pole contribution is changed to the off-shell power contribution, while T eq is the one 
when the two temperature dependent terms in the stationary number density 
compete equally; 

Importance of the off-shell effect is determined by the remnant density at T = T eq ; 
(rz/T 3 ) T=Tcq . Besides a weak logarithmic dependence, this quantity has the factor, 
T/M. The effect is thus more pronounced for a larger rate T, or a stronger coupling. 



4. Some application and extention 

Physical processes that occur after this transient time are prolonged, compared 
to the instantaneous decay in the exponential law. One of the important effects of 
this sort is baryogenesis at the GUT epoch. In the scenario of a heavy X boson 
decay the out-of-equilibrium condition was considered to impose a rather severe 
constraint on particle physics and cosmology. Q This constraint arises due to that 
the abundance of X and X bosons should not be too much suppressed when the 
asymmetry generation starts. The on-shell kinematics requires that the temperature 
Ta at the Hubble H = F must not be too much larger than the mass of X boson; 



10 



thus T d < m x . Otherwise, the inverse decay can take place frequently, to suppress 
the number density of X boson. This leads to the lower bound of the X boson mass, 
m x > 0[1] x « 10 16 GeV. 

This constraint is however based on the on-shell kinematics; thermal quarks 
and leptons are energetically difficult to produce heavy X bosons when the cosmic 
temperature T < mx- But, as we discussed here, heavy X bosons that were once 
abundant do not disappear immediately due to the off-shell effect. Thus, one can 
anticipate that the above lower bound of the X mass is considerably relaxed. |TJ For 
details of this problem, I refer to our forthcoming paper, [g] 

As an extention of the model so far discussed, it is instructive to treat the case 
of pair annihilation model. The model can describe two particle annihilation into 
lighter particles, 

ci c 2 -> b(u) , 

again &t(u;)|0) taken to be two-body states of light particles. Its model Hamiltonian 
is 

H = uji c\c\ + uj 2 c\c2 + / dujuj b\uj)b{uj) + / duj\Ja[uj) ( b\u)ciC2 + (h.c.) ) . 

Jlu c J UJ C 

(34) 

This time the model does not allow exact solution, but one can use the method 
of Feynman and Vernon ||, to integrate the effect of environment variables. We 
leave details of this technical part to our future paper, but let me discuss how one 
may derive a Boltzmann-like equation. The key for this is to neglect a long-time 
correlation and to take the t — > oo limit in order to erase the off-shell effect. After 
cutting off the long-time correlation, we find that 

-(cl(t)ci(t)) = / duja(uj) ds 



dt Juj c Jo 

2cos(lu-uj 1 -uj 2 )s ( - /i/ 2 (/th + 1) + (A + l)(/2 + l)/th ) ) , (35) 

where /th(^) = 1/ ( e ^ w ~ 1) an d fi{oJi) are occupation numbers of i. One may then 

lim / ds 2 cos(ciJ — oj\ — uj 2 )s = 2tt 5{uj — uj% — cj 2 ) , (36) 
;^oo Jq 



use 

rt 



t— >oo 

to get 



j t (4(tMt)} = r ( - /iM/aMC/th + 1) + (/i(^i) + l)(/ 2 (wa) + l)/ th ) , (37) 



If 



with f th = /th(^ = W\ +0J2) ■ This is the Boltzmann equation with the rate averaged 
in thermal medium. Our formalism can thus describe the kinetic equation for the 
pair annihilation in thermal medium, once the thermal distribution of two-body 
states /th(^) is given. 
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